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AN XMM-NEWTON STUDY OF THE 401 HZ ACCRETING PULSAR SAX J1808.4-3658 IN QUIESCENCE 
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ABSTRACT 

SAX J1808. 4-3658 is a unique source being the first Low Mass X-ray Binary showing coherent pulsations at 
a spin period comparable to that of millisecond radio pulsars. Here we present an XMM-Newton observation of 
SAX J1808. 4-3658 in quiescence, the first which assessed its quiescent luminosity and spectrum with good signal 
to noise. XMM-Newton did not reveal other sources in the vicinity of SAX J1808.4-3658 likely indicating that 
the source was also detected by previous BeppoSAX and ASCA observations, even if with large positional and 
flux uncertainties. We derive a 0.5-10 keV unabsorbed luminosity of Lx = 5 x 10 31 erg s~\ a relatively low 
value compared with other neutron star soft X-ray transient sources. At variance with other soft X-ray transients, 
the quiescent spectrum of SAX J1808.4-3658 was dominated by a hard (T ~ 1.5) power law with only a minor 
contribution ( ^ 10%) from a soft black body component. If the power law originates in the shock between the 
wind of a turned-on radio pulsar and matter outflowing from the companion, then a spin-down to X-ray luminosity 
conversion efficiency of 77 ~ 10~ 3 is derived; this is in line with the value estimated from the eclipsing radio pulsar 
PSR J1740-5340. Within the deep crustal heating model, the faintness of the blackbody-like component indicates 
that SAX J1808.4-3658 likely hosts a massive neutron star (M > 1.7 M ). 

Subject headings: accretion, accretion disks - 
neutron 



binaries: close — star: individual (SAX J1808-3658) — stars: 



1. INTRODUCTION 

Neutron star Soft X-ray Transients (SXRTs) when in out- 
burst closely resemble persistent Low Mass X-ray Binaries 
(LMXRBs). In the last few years it has become clear that 
SXRT sources form a rather inhomogeneous class (see Cam- 
pana et al. 1998a for a review) comprising sources with well 
defined outbursts as well as sources with long on/off activity 
periods. Moreover, sources displaying bright outbursts with 
peak X-ray luminosities Lx ~ 10 37 — 10 38 erg s _1 appear 
to be different from sources showing only faint outbursts reach- 
ing Lx ~ 10 36 — 10 37 erg s _1 , especially in the Galactic center 
region (Heise et al. 1998; King 2000; in't Zand 2001a). 

A major leap forward came with the discovery of SAX 
J1808.4-3658, a bursting SXRT reaching a maximum luminos- 
ity of ~ 2 x 10 36 erg s _1 (for a distance of 2.5 kpc; in'tZand al. 
2001b). In April 1998 the source resumed activity and RXTE 
observations revealed coherent ~ 401 Hz pulsations, the first 
detected in the persistent emission of a neutron star LMXRB. 
These testify to the presence of magnetic polar cap accretion 
onto a fast rotating magnetic neutron star (Wijnands & van der 
Klis 1998; Chakrabarty & Morgan 1998). The inferred mag- 
netic field strength of SAX J1808.4-3658 is in the 10 8 - 10 9 G 
range (see Psaltis & Chakrabarty 1999), providing convincing 
evidence for the long suspected LMXRB -millisecond pulsar 
connection. 

SXRTs spend much of their time in quiescence. The origin of 
the quiescent X-ray emission is still uncertain. In the last few 
years several sources have been studied in detail; the picture 



that emerged is that the quiescent spectrum consists of a soft 
component plus a high energy excess. The former component 
is often fit with a blackbody spectrum with an equivalent radius 
of 1-2 km and temperatures in the 0.1-0.3 keV range or with a 
neutron star atmosphere model with an equivalent radius con- 
sistent with the entire neutron star surface and slightly smaller 
temperatures (Brown et al. 1998; Rutledge et al. 2000). The 
high energy component is well represented by a power law with 
photon index r ~ 1 — 2 (Asai et al. 1996, 1998; Campana et 
al. 1998b, 2000). In all sources observed so far the quiescent 
luminosity ranges between 10 32 — 10 33 erg s _1 , indicating a 
clear difference with black hole transients in quiescence that 
have a lower X-ray luminosity (Garcia et al. 2001 ; Campana & 
Stella 2000). Flux and spectral variability have been reported in 
Aql X-l and KS 1731-260 during quiescence (Campana et al. 
1997; Rutledge et al. 2002; Wijnands et al. 2002a). This poses 
severe limitations on the emission mechanisms responsible for 
the quiescent luminosity. 

Among neutron star SXRTs, SAX J1808.4-3658 stands out 
for having, while in outburst, a magnetosphere and, of course, 
an accurately measured spin period. These characteristics make 
the source very well suited for testing the predictions of models 
for the quiescent emission in which the presence of a sizeable 
magnetic field plays a crucial role (Stella et al. 1994). SAX 
J1808.4-3658 was detected in quiescence (Stella et al. 2000; 
Dotani et al. 2001; Wijnands et al. 2001) though with large un- 
certainties. Here we report on an XMM-Newton observation of 
SAX J1808.4-3658 in quiescence, the first to detect the source 
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with a good signal to noise ratio. 

2. DATA ANALYSIS 

SAX J1808.4-3658 was observed on March 24 2001 with 
XMM-Newton EPIC, consisting of two metal oxide semicon- 
ductor (MOS) cameras (Watson et al. 2000) and one pn camera 
(Struder et al. 2000). Medium filters were used for the MOS 
cameras and the thin filter for the pn camera. The pn cam- 
era was operated in timing mode in order to search for pulsa- 
tions in case the source were sufficiently bright. SAX J 1808.4- 
3658 turned out to be faint (see below), preventing meaningful 
searches for pulsed emission. We do not discuss the pn camera 
data in the following and concentrate on data from MOS cam- 
eras that were operated in full frame mode with a read-out time 
of 2.6 s. 

We extracted the event file starting from the raw data using 
the Standard Analysis System (SAS) version 5.3.0. Data were 
manually screened to remove any remaining bright pixels or hot 
columns. Periods in which the background was high because of 
soft proton flares were excluded using an intensity filter: we re- 
jected all events accumulated in the external CCDs (2-7) when 
the total count rate exceeded 15 counts in 100 s in the 10-12.4 
keV band for the MOS1 and MOS2, independently. We ob- 
tained a net exposure time of 33.2 ks and 34.4 ks for MOS1 
and MOS2, respectively (this is slightly different from what re- 
ported in Campana et al. 2002, due to the use of standard pro- 
cessed files with SAS 5.0.3). Spectra were accumulated sepa- 
rately for MOS1 and MOS2. Event grades higher than 12 were 
filtered out. Photons were grouped within XMMS ELECT to the 
nominal MOS resolution of 15 eV. 

2.1. Source position 

Previous observations aiming at determining the position of 
SAX J1808.4-3658 were performed with BeppoSAX (Stella et 
al. 2000; Wijnands et al. 2002b) and ASCA (Dotani, Asai 
& Wijnands 2000; see Table 1), but due to the limited angular 
resolution of these satellites the resulting picture was quite con- 
fused. In particular it was unclear whether the source revealed 
at a low flux of ~ 10~ 13 erg s _1 cm~ 2 and with large posi- 
tional uncertainties was indeed SAX J1808.4-3658 (see Fig. 1 
and Tab. 1). In fact, Wijnands et al. (2002b) claimed that SAX 
J1808.4-3658 was only detected in the ASCA data of Septem- 
ber 1999, while neither BeppoSAX observations detected the 
millisecond pulsar but rather revealed a nearby source, SAX 
J1808.6-3658. 

The deep image we obtained with EPIC allowed us to clarify 
this confused picture by revealing several faint sources in the 
region. Source detection was obtained following the prescrip- 
tions of Baldi et al. (2002). This revealed 20 sources in the 
central part (12' x 9') of the field of view. The brightest one 
was clearly coincident with the radio position of SAX J 1808.6- 
3658. A few other sources were contained in the error regions 
of sources previously seen with ASCA and in the longer Bep- 
poSAX observation (see Fig. 1). 

A faint source is visible about 2' East of SAX J1808.4-3658 
(RA(J2000): 18 h 08 m 36 s .6; DEC(J2000):-36° 58' 01"). This 
is within the error region of the source detected by Wijnands 
et al. (2002b) in the March 2000 BeppoSAX observation. If it 
had the same luminosity that we measure with EPIC (a count 
rate a factor ~ 5 lower than that of SAX J1808.4-3658), it was 
probably too faint to be detected by BeppoSAX, but a small 
contamination from this source could perhaps explain why all 



the error regions for SAX J1808.4-3658 derived with low an- 
gular resolution instruments are systematically centered to the 
East of the radio position. A further bias is provided by a 
2.9' to the brighter source East (RA(J2000): 18 h 08 m 42 s .2; 
DEC(J2000):-36° 58' 44") which has a comparable flux to 
SAX J1808. 4-3658 and went undetected in previous observa- 
tions (Wijnands et al. 2002b). Furthermore, given the source 
density revealed in this field by XMM-Newton, the probability 
of finding by chance at least one source within the BeppoSAX 
90% error circle of SAX J1808.6-3658 is - 70%. It appears 
to be unlikely that this faint source was much brighter and at 
the same time SAX J1808.4-3658 much fainter just during the 
BeppoSAX observations and conclude that SAX J1808.4-3658 
was the source detected both in March 1999 and March 2000. 

2.2. X-ray spectrum and luminosity 

In total we collected 91 and 77 counts with the MOS1 and 
MOS2 detectors, respectively, within a 10" radius circle cen- 
tered on the position of SAX J1808.4-3658. This corresponds 
to about 55 — 60% of the source flux for energies between 1- 
5 keV (extending the extraction region would just decrease the 
signal to noise ratio, due to the background and source faint- 
ness). The background was extracted from an annular region 
with inner and outer radii of 0.5' and 2.5' (excluding contam- 
inating sources), respectively. Within the SAX J1808.4-3658 
extraction region, about 15 — 20% of the counts are to be at- 
tributed to the local background. 

The MOS 1 and MOS2 spectra were rebinned to contain at 
least 15 photons per channel. The spectral analysis was carried 
out with XSPEC v.l 1.1.0. We adopted the latest on-axis stan- 
dard response matrices (ml_medv9q20t5r6_all_15 . rsp 
and m2_medv9q20t5r6_all_15 . rsp) and correct the 
fluxes for the fraction of photons collected within the extrac- 
tion radius. The spectral analysis was carried out in the 0.3-7 
keV energy range (Kirsch private communication and 2002). 
Given that the cross calibration between MOS1 and MOS2 
agrees within 5% (Kirsch 2002) we did not include a con- 
stant for the two MOS instruments. The statistics was rela- 
tively poor and the spectrum could be well fit by an absorbed 
power law (or a bremsstrahlung) model (see Table 2 and Fig. 
2). A power law model gave a better fit (x 2 od = 1-1) w i tn 
r = 1.51q'3 (90% confidence level) and a column density of 
N H = 0.3 x 10 21 cm- 2 (< 1.1 x 10 21 cm- 2 90% c.L, we used 
here the Wisconsin cross sections in the wabs model). We note 
that a model with the column density fixed at the value derived 
from the source data in outburst (e.g. Gilfanov et al. 1998) and 
consistent with the galactic value of 1.3 x 10 21 cm -2 , provides 
also a good fit with T = 1.8 ± 0.3 and Xr cd = 1-2 (with a null 
hypothesis probability of 20%). The unabsorbed 0.5-10 keV 
luminosity of SAX J1808.4-3658 was 5 x 10 31 erg s" 1 for the 
power law model and a (revised) distance of 2.5 kpc (in't Zand 
et al. 2001b). Single component black body or neutron star at- 
mosphere models (Gaensicke, Braje & Romani 2002) did not 
provide an adequate fit to the data with x 2 od > 3. 

Given that the best studied quiescent SXRTs display a soft 
black body component plus a hard power law tail we con- 
sider also this model (even though strictly not required by 
the data). In this case we fixed again the column density to 
1.3 x 10 21 cm- 2 . The power law fit alone could adequately 
describe the data so that we first fixed the best fit power law 
and derived constraints on the black body component. SXRTs 
in quiescence display black body temperatures ranging between 
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0.1-0.3 keV (e.g. Rutledge et al. 2000). By using this range 
we derived a 90% upper limit to the normalization of the black 
body component and in turn on the (unabsorbed, bolometric) 
soft black body flux of ~ 2 x 10~ 15 erg s _1 cm~ 2 . This cor- 
responds to a factor of 15 less than the power law luminosity 
in the 0.5-10 keV band. If we repeat the same exercise with 
a neutron star atmosphere model (Gaensicke, Braje & Romani 
2002) we derive a factor of 30 less luminosity in the soft com- 
ponent. Alternatively, by fixing the normalisation ratio of the 
black body and power law components as observed in the best 
studied quiescent SXRTs (e.g. Campana et al. 1998b, even if 
recent observations on Aql X-l challange this picture, Rutledge 
et al. 2002) we obtained a very low black body temperature. By 
fixing also the column density to 1.3 x 10 21 cm -2 we succeeded 
in obtaining a good fit x 2 od = 1.2 with a black body tempera- 
ture of kT = 0.161^5 keV and a power law V = 1.5 ± 0.5. 
The source flux is in any case low 2.5 x 10 31 erg s^ 1 . With 
a neutron star atosphere model we obtain an equally good fit 
( x 2 ed = 1.2) with a temperature k T = 0.09±g;2| keV and an 
emitting radius R < 1.1 km. 

3. DISCUSSION 

The high throughput and good angular resolution of XMM- 
Newton provided the first firm determination of the quiescent 
luminosity of SAX J1808.4-3658. This is lower by a factor of 
~ 2 than previous best fit estimates. However, once the large 
uncertainties in the spectral parameters are taken into account, 
a fairly constant luminosity is inferred since March 1999 at a 
level of 5 x 10 31 erg s _1 (at a distance of 2.5 kpc; note that the 
reduced \ 2 of a fit with a constant is 1.8, 15% null hypothesis 
variability, see Fig. 3). We conclude that the source variability 
issue raised by Dotani et al. (2000) is questionable. The source 
luminosity is a factor of two lower than that usually measured in 
SXRTs in quiescence; this ranges between 10 32 — 10 33 erg s _1 
(e.g. Campana et al. 1998a; Campana & Stella 2000; Garcia et 
al. 2001). For a distance of 4 kpc, however, it would be fully 
consistent with what usually observed in other SXRTs. 

More interestingly, if we consider the spectral fit with a 
power law component plus a black body, the soft X-ray com- 
ponent comprises only a small part ( ;$ 10%) of the total lu- 
minosity, whereas in the great majority of SXRTs it accounts 
for about half of the luminosity in the 0.5-10 keV energy band 
(e.g. Campana et al. 1998b, 2000; see however Rutledge et al. 
2002). The soft component is usually ascribed to the cooling 
of the neutron star surface powered by the deep nuclear heat- 
ing that the neutron star receives during each outburst (Brown 
et al. 1998; Rutledge et al. 2000; Colpi et al. 2001; see also 
Campana et al. 1998a). RXTE All Sky Monitor (ASM) pro- 
vides us with a continuous monitoring of the high level activity 
of SAX J1808. 4-3658. From these data we can extrapolate a 



mean mass transfer rate of ~ 5 x 10~ 12 M yr _1 , in line with 
previous estimates (Bildsten & Chakrabarty 2001). This mean 
mass inflow rate translates into a soft quiescent luminosity of 

10 32 erg s _1 within the 0.5-10 keV energy band from deep 
crustal heating (Brown et al. 1998; Colpi et al. 2001). This 
is a factor of ~ 10 higher than observed. Given the fact that 
predictions from deep nuclear heating are quite robust, one is 
led to conclude that an additional source of cooling is present. 
A simple and well known solution is when the direct Urea pro- 
cess is allowed in the neutron star core and, in turn, neutrino 
cooling does affect the neutron star thermal evolution (Colpi et 
al. 2001). This can occur only for massive neutron stars with 
masses higher than 1.7 — 1.8 M Q . If this interpretation were 
correct the neutron star of SAX J1808. 4-3658 have to be fairly 
massive, in agreement with accretion spin-up scenarios. 

The main contribution to the quiescent luminosity of SAX 
J1808. 4-3658 appears to derive from the power law compo- 
nent. A pure propeller contribution is ruled out since this 
mechanism should stop operating at a luminosity of about 

10 33 erg s _1 , with the turing on of a radio pulsar (e.g. Cam- 
pana et al. 1998a, 1998b). One interpretation for this power law 
relies on the emission at the shock front between the relativistic 
wind of a radio pulsar and matter outflowing from the compan- 
ion star (see the discussion in Stella et al. 2000). Being the 
first SXRT for which the presence of a sizable magnetic field is 
unambiguously established, SAX J1808.4-3658 can be used to 
infer the efficiency rj with which spin-down power is converted 
into 0.5-10 keV luminosity. We derive r\ <~ 5 x 10~ 3 B\ , where 
B s = B/10 8 G is neutron star magnetic field (as derived from 
the accretion luminosity at the propeller onset, Gilfanov et al. 
1998; see also Psaltis & Chakrabarty 1999). 

Recently, an ideal laboratory for studying the shock emis- 
sion has been discovered: this is PSR J1740-5340, a 3.7 ms 
radio pulsar in a 32.5 hr orbit around a Roche lobe filling main 
sequence companion, which display partial and total eclipsed 
over a wide range of orbital phases (D'Amico et al. 2001; Fer- 
raro et al. 2001; Burderi, D'Antona & Burgay 2002). The 
X-ray luminosity of PSR J1740-5340 measured by Chandra 
(8 x 10 30 erg s~\ 0.5-2.5 keV range, unabsorbed; Grindlay 
et al. 2001) implies r\ ~ 10~ 4 . Taking this system as a proto- 
type for modelling the shock emission power (which is likely 
less efficient since the radio pulsar is not completely engulfed) 
and scaling the efficiency as the square of the orbital separation, 
one infers f7^3xl0~ 3 , well consistent with the value above, 
given the uncertainties involved. 

We thank A. Baldi for the source detection algorithm and A. 
De Luca for the image hot pixels cleaning procedure. We ac- 
knowledge useful comments by an anonymous referee which 
helped improving the paper. 
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Table 1 
SAX J1808. 4-3658 positions. 



Obs. date 


Satell. 


Exp. time 


Error radius 


Distance* 


Refs.T 


Mar. 1999+ 


BSAX 


20 ks 


1.50' (90%) 


1.43' 


1 


Sep. 1999 


ASCA 


61 ks 


0.24' (68%) 


0.19' 


2 


Mar. 2000 


BSAX 


89 ks 


0.55' (68%) 


1.63' 


3 


Mar. 2001 


XMM 


33 ks 


0.10' (90%) 


0.04' 


4 



* Distance of the detected source from the radio position of SAX J1808.4-3658. 

t 1: Stella et al. (2000); 2: Dotani et al. (2000); 3: Wijnands et al. (2002); 4: this paper. 

+ Source position derived is affected by a large uncertainty due to source faintness and unfavorable roll angle resulting in a systematic 
uncertainty of about 35" on top of the intrinsic positional error (Perri & Capalbi 2002). 



Table 2 

Spectral models and related luminosities. 



Model 


N H 
(10 21 cm" 2 ) 


TlkT 


Xrcd 


Flux 

(10- 14 cgs) 


L x 

(10 31 erg s" 1 ) 


Black body+ 


OQ+0.3 


0.52i^ keV 


3.2 


1.5 


2.2 


Power law 


0.3+ ' 7 


1 5 +0 - 2 


1.1 


2.9 


4.6 


Bremsstrahlung 


00 +o.9 


15.5_io.4 keV 


1.1 


2.8 


4.4 


NS atmosphere* 


0.4+ 11 


0.27±^° 2 keV 


4.4 


1.1 


1.8 



Values without an explicit error indicate that the parameter is unconstrained. Fluxes and luminosities are in the 0.5-10 keV energy 
band. Fluxes are absorbed and luminosities are unabsorbed. 

+ Even fixing the column density to 1.3 x 10 21 cm~ 2 we derive a temperature 0.5 keV and a reduced x 2 = 3.2. 

* Model for a hydrogen atmosphere by Gaensicke, Braje & Romani (2002). The temperature collapses to the maximum allowed 

value of 0.275 keV. 
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FIG. 1 . — Upper panel: XMM-Newton observation of the SAX Jl 808.4-3658 field with overlaid sources detected during previous observations with BeppoSAX 
and ASCA (see text for more details and Table 1). Lower panel: Same field as above with XMM-Newton detections and SAX J1808. 4-3658 previous detections. 
The error circles for the XMM-Newton sources have been enlarged to 15" for display purposes. 
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FIG. 2. — MOS1 (open dots) and MOS2 (filled dots) spectrum of SAX J1808.4-3658. Overlaid is the fit with an absorbed power law model described in the text. 
In the lower panel are reported the residuals of the fit. 
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FIG. 3. — Quiescent 0.5-10 keV unabsorbed luminosities of SAX J1808.4-3658 observed with BeppoSAX, ASCA and XMM-Newton with relative errors (1 a). 
The dashed line represents the fit with a constant value. In the case of the BSAX2 observation a smaller uncertainty was reported by Wijnands et al. (2002b). The 
small number of counts, the high BeppoSAX background within the extraction radius (comparable to the source flux), the a posteriori discovery of a weak source 
within this radius make the flux uncertainty quoted by Wijnands et al. (2002b) too small. 



